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ABSTRACT 

In this paper we prove the existence of one-sided filters, for spectral Fourier approxima- 
tions of discontinuous functions, which can recover spectral accuracy up to the discontinuity 
from one side. We also use a least square procedure to construct such a filter and test it on 
several discontinuous functions numerically. 
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1 Introduction 


A spectral Fourier approximation for a 2i r periodic function u(x) is 

N 

u n(x) = U k e ikx (1.1) 

k=-N 

where the Fourier coefficients u k are defined by 


1 r 2ir 

u k = — I u(x)e~ ,kx dx 
2 tt J o 

for Fourier Galerkin, and by 


( 1 . 2 ) 


u k = 


1 

2N+ 1 


2N 

E u (xj)e- ik *i, 


i= o 


x,= 

’ 2N+1 


(1.3) 


for Fourier collocation. Spectral methods are well-known for their high accuracy in the ap- 
proximation of smooth functions and in solving partial differential equations with smooth 
solutions [3], [2]. For a discontinuous function, however, unmodulated spectral approxima- 
tion produces Gibbs oscillations and yields first order accuracy even in the smooth region 
[3]. The good news is that, even if the accuracy is poor in the point-wise sense, it is still 
excellent for the moments [4]: 

Lemma 1.1 If u(x) £ L 2 [0,27r] and ujv(x) is its Fourier Galerkin sum (1.1)-(1.2), then 
for any 27r periodic function v{x ) £ C°°, we have 


f 2 * 

I (u(x) — ujv(x))t)(x)da 

J o 


< 


Ill’ll 




N > 


|u||x,j Vs > 1 


(1.4) 


Proof : By the orthogonality of trigonometric polynomials, we have 


y2ir 

I (u(x) — Uff(x))v(x)dx 

J 0 


< 


I J o {u(x) - u N (x))(v(x) - va r(x))dx 
||m||l» ||v — ^Ar|jx, 2 


Integration by parts yields 


1 



v k = — r v(x)e~ ikx dx = -V ■ r v^(x)e~ ikx dx 
27r Jo k j 2t r(-zJfc)* Jo v ' 


We then have 


|u - vat ||z,a 



|v w ||l 3 

N' 


□ 


Lemma 1.1 implies that the spectral approximation u^(x) does indeed contain the in- 
formation of u(x ) by preserving accurately the moments of u(x) against any C°° functions. 
Various filters axe designed in the literature to extract this information. The early work in 
this direction includes [8] by Majda, McDonough and Osher and [7] by Kreiss and Oliger. 
The filters, in the phase space, are of the form 

u n ( x ) = X) < 7 f«* e ‘ fce (1*5) 

k=-N 

where < 7 k = a^ k are real numbers, which decay smoothly from 1 to 0 when |&| goes from 0 
to N. The filter (1.5) is also equivalent to a convolution in the physical space 

1 f 2 * 

u-n( x ) = u N * K N (x) = — / u N (y)K N (x - y)dy (1.6) 

Z7T JO 

with a kernal Kn{x) defined by 

Kn(x) = £ (1.7) 

k~-N 

Vandeven [10] studied this type of filters in more detail. He proved the following result, 
which we will use in Section 2. For simplicity of presentations and without loss of generality, 
we assume the function u(x) has only one discontinuity located at x = 0, i.e., u(x) is smooth 
but not periodic over [0,27r). 
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Lemma 1.2 (Vandeven) If u(x) is an analytic but not periodic function in [0,27 t), then 
for any 0 < £ < 1, the filtered Fourier sum u^(x) defined by (1.5) with 


p=Ni 


satisfies the following error estimates: 


( 1 . 8 ) 


} max |u(x) — u^-(x)| < 


N e 




(CN*) Ni 


where C and /? are constants independent of N. 


(1.9) 


Similar results were also obtained in Gottlieb and Tadmor [4], by directly constructing 
the filter kernel K^{x ) in (1.6). 

In this paper, we will use C or C for a generic constant independent of N, which may be 
different in different locations. 

Lemma 1.2 establishes the spectral convergence of the filtered Fourier sum (1.5) to the 
function u(x), in a distance jynrr away from the discontinuity x = 0 ( mod 2 tt). If u(x) 
is not piecewise analytic but piecewise C T , similar estimates can be obtained for algebraic 
convergence. 

The choice of in (1.8) is not unique. In practice, exponential filter of the form 

aS = ( 1 . 10 ) 

is often used due to its simplicity and good numerical results. Here 2p is the order of the 
filter and can be chosen depending on N. a is a constant chosen so that e~ a , the filter for 
the last mode, is machine zero. Other frequently used filters include the raised cosine filter 
and the Gottlieb-Tadmor filter [4]. We refer the readers to Kopriva [6] for the definition and 
an extensive comparison of various filters in practical calculations. 
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All of the filters mentioned above are two-sided filters. That is, a two-sided region around 
the discontinuity has to be excluded in the error estimates. See for example (1.9). This is 
expected since one uses a symmetric kernel Kn(x) in (1.6), as a result of using real numbers 
a k = a^ k in (1.5). The Gibbs oscillation is still present in the filtered sum u < ^(x) near the 
discontinuity, see Figure 1. This is usually unacceptable for solving nonlinear partial differ- 
ential equations since these oscillations may eventually pollute the smooth regions and/or 
trigger nonlinear instability. In [1], we proposed a way to overcome this difficulty by intro- 
ducing non-smooth functions (saw-tooth functions) into the basis of trigonometric functions. 
It would be nice, however, if one could work within the trigonometric polynomial basis and 
use filters to completely remove the Gibbs oscillations and to obtain uniform convergence up 
to the discontinuity. This is the motivation for the one-sided filters discussed in this paper. 
We refer the readers to Mock and Lax [9] for the early work in this direction. 

In Section 2 we prove the existence of one-sided filters, i.e., we prove the existence of 
complex numbers c k , such that for any analytic but not periodic function u(a:) in [0, 27r), 
the filtered Fourier sum (1.5) satisfies a uniform error estimate 


where 


max |u(x) — u%t(x ) I < 

0<x<iii v ' r<\ — 


N 0 

(CNi)"* 


( 1 . 11 ) 


xr = 2tt - 2 , (1-12) 

JV 1_ 4* 

and C and (3 axe again constants independent of N . This filter is naturally labelled “right- 
sided” filters since it can recover the spectral accuracy up to the discontinuity x = 0 from 
the right. A symmetry consideration leads to the “left-sided” filters by taking the conjugate 
of the complex numbers a k . 

The proof in Section 2 is constructive. However, the filters obtained that way can not 
be satisfactorily applied in actual computations unless N is extremely large. For a practical 
range of N between 8 and 32 (16 to 64 grid points for collocation), we try to find a better 
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one-sided filter through a least-square procedure described in Section 3. Several numerical 
examples are also provided in that section. 

We point out that the one-sided filters discussed in this paper are neither unique nor 
necessarily optimal in practical calculations. Work is under way to explore the construction 
of one-sided filters directly in the physical space (1.6). A systemetic investigation about 
one-sided filters in the phase space, similar to the approach used in [10], is also under 
consideration. 


2 Existence of One-Sided Filters 


The main result in this section is the following theorem: 

Theorem 2.1 If u(x) is an analytic (but not periodic) function in [0, 27 t), then for any 

0 < e < y, the one-sided filter defined by 


where 




( 2 . 1 ) 


A = 


1 

N 1 -*’ 


771 = N* 


( 2 . 2 ) 


and is the Vandeven two-sided filter defined by (1.8), produces a filtered Fourier sum 


(1.5) satisfying 


N p 

where xr is given by (1.12), C and /3 are constants independent of N . 
Proof : We have, by (2.1) and (1.5), 


(2.3) 


Utf(x) 


E <^e( 7 ) (-1) ,+1 ^ A ^ 

k=-N 1=1 V / 
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(2.4) 


= E ( ? ) (-1) ,+, “» T (- + JA) 

= E(7)(- 1 )' +1 “( I+iA ) + £ ‘ 

where Ei = ( 7 ) (~l) ,+1 0*Jf r (* + ^A) — u(x + /A)), and 


u N T i x )= E °l' N Uke' kx (2.5) 

k=-N 

is the two-sided filtered Fourier sum. By (2.2) and Lemma 1.2, E\ is small whenever 0 < 
x < xjj, i.e., whenever < x + IA < 2n — for 1 < l < m: 


is, I < E 


m 


N 0 


= 2 " 


N 0 


N 0 


<5 = 2 


SV 1 1 (CNiy* (CJV5)W* (<7atJ)»« ’ ~ 2 

On the other hand, by Newton’s formula, the sum on the last line of (2.4) satisfies, for 
some x < £ < x + mA, 


E [ i ) (-i) I+ v* + w = u( X ) - »w(o(-a r 

Since we assume u(x) is analytic in [0, 27r), we have, for some p > 0, 

|u< m >(x)| < M(p)- 


,m! 


We then have, by (2.2) and the Stirling’s formula lin^^oo m "^(£Z rn = 1, 


Cm! 


N» 


Nt 


| tt (-)(f)A m | < — __l_ — < 111 < 

K } - (pW 1 -*)" 1 - (CiV 1 -?*)"* “ (CNiyrf 

where the last inequality is valid because £ < This completes the proof. 


Remark (1) The filter (2.1) is based on the following simple idea: since the two-sided filtered 
Fourier sum u^ T (x) in (2.5) is a good approximation to u(x) in a region A away from the 
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discontinuity, we can use extrapolations from points inside this region to approximate points 
outside of it. The right-sided filtered Fourier sum ujr(x) based on (2.1) is just such an 
extrapolation using u^ T (x). In practice, this simple, equally spaced extrapolation may not 
be the best. One would expect Chebyshev extrapolations to do a better job. 

(2) A left-sided filter is obtained by taking the complex conjugate of a k in (2.1). The 
proof follows the same lines. 

3 Efficient One-Sided Filters and Numerical Exam- 
ples 

The one-sided filters described in Section 2 are good asymptotically, but a very large N is 
needed in order for m in (2.2) to be of reasonable size. It is our experience that asymptotically 
correct choices are often not necessarily optimal for small N . We thus use a least square 
procedure, described below, with the objective of obtaining more efficient one-sided filters 
for a practical range of N between 8 and 32 (between 16 and 64 grid points for collocation): 


Least Square Procedure: We make an ansaze 


„N _ . 

a k = < 


f + *>o 

* < 0 


(3.1) 


T N 

r -Jt> 


where the weight function w{t) is defined by 


w(t) = e-^tby (3.2) 

al ,N is the Vandeven two-sided filter defined by (1.8) with a parameter p, and q(t) € P ,, the 
collection of all polynomials of degree at most s, with complex coefficients, defined on [0, 1]. 
The right-sided filter o k is then defined by taking the following minimum: 


min 

q(t)eP, 


in £ /V)W*) - *T<k 

eP. ^ Jo 


(3.3) 
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where A is a parameter between 0 and 27 r, and (x r )^(a:) is the filtered Fourier sum (1.5), 
using the filter (3.1), of the function u(x ) = x T . 

□ 

The parameters involved in the Least Square Procedure are p for the order of accuracy 
of the two-sided filter (1.8), s for the degree of polynomials q(t) in (3.1), R for the number of 
terms in the summation (3.3), and A in (3.3) for the domain of integration. The implemen- 
tation of the minimization (3.3) involves the choice of a basis for P s (we use a Chebyshev 
basis over [0, 1]) and the solution of the resulting linear system. 

The condition a^ k = a k is imposed in (3.1) in order to get a real kernel in (1.7) or, 
equivalently, to get a real filtered Fourier sum (1.5) for a real function u(x). The exponential 
weight function (3.2) is used so that the term added to the two-sided filter <jJ ,n will not 
destroy the accuracy (see [10] for details of the relation between properties of c k and the 
accuracy of the filters). The minimization (3.3) is based upon the assumption that the right- 
sided filtered Fourier sums for u(x) = x r , with 1 < r < R, should be uniformly accurate in 
the interval [0,A]. We have also tested the procedure by replacing x r in (3.3) by the r-th 
order Chebyshev polynomial T r (x), obtaining similar results. 

A crucial issue for the success of this approach is the sensitivity of the filters obtained 
with respect to the paramters p, s, R and A. According to the experience for two-sided 
filters, p should be chosen related to N. One could also try the (very costly) procedure of 
minimization over some of those paramters in certain ranges. The filter might be expected 
to work well for polynomials because of the choice u(x) = x T in the minimization (3.3) but 
must be extensively tested for other functions. 

We perform numerical experiments using different parameters p, s, R, A to get the filters, 
then testing them on a non-polynomial function 

u(x) = cos(1.5x) iE[ 0, 27 t) (3-4) 
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Since we enforce 2ir periodicity, this function and all its derivatives are discontinuous at 
x = 0 ( mod 2 tt). We have experimented with p between 3 and 15 depending on the size 
of N (the choice is based on our experience with two-sided filters). It turns out that the 
result is most sensative to the choice of s and the optimal s increases rather rapidly with 
N. This indicates that the choice for the polynomial space P, or the exponential weight 
function in (3.2) may not be adequate for large N . It also limits the size of N we can use in 
this approach since the resulting linear system, which we solve by using a routine in IMSL, 
becomes ill-conditioned for large s. However, the result seems not very sensative to the 
choice of R (we have tried 4 < R < 15) and A (bounded away from 0 and 27 r). 

We now show the results obtained for N = 8 with 3 = 2, p = 3, R = 5 and A = 7r; for 
N = 16 with 3 = 5, p = 5, R = 5 and A = 7r; and for N = 32 with s = 18, p = 6, R = 6 and 
A = 7r, in Figures 2 to 6. 

In Figure 2, we plot the filter aj? for its real and imaginary parts. No special pattern can 
be observed. In Figure 3, we plot the filter for its magnitudes and arguments. Clearly it 
shows a straight line for the arguments immediately after the first few modes. This interesting 
phenomena can be loosely explained in the following way: A straight line for the arguments 
corresponds to a pure shifting for the kernel (1.7). It makes the kernel non-symmetric around 
zero and mainly supported in one-side, allowing for one-sided recovery of accuracies. For the 
first few modes the arguments have to be close to zero to ensure accuracy of the filter for 
the low modes. 

In Figure 4, we plot the filter kernels (1.7) in the physical space. We can see that they 
are indeed one-sided approximate 8 functions, i.e., they are supported mainly to the left side 
of the discontinuity 0 ( mod 2i r). 

In Figure 5, we plot the one-sided filtered Fourier sum (1.5), using plus signs, against 
the exact solution (3.4) (shifted by 7r to show more clearly the discontinuity), with N = 16. 
Each point is either left-side filtered or right-side filtered according to which side it is closer 
to the discontinuity, we can clearly see the advantage of using one-sided filters in getting 


9 


accurate, non-oscillatory results. 

Finally, in Figure 6, we plot the errors of the one-sided filtered Fourier sum (1.5) from 
the exact function (3.4), with N = 8, 16,32, on a logrithm scale. Again, each point is either 
left-side filtered or right-side filtered according to which side it is closer to the discontinuity. 
We can see the uniform convergence up to the discontinuity and a faster than algebraic 
convergence for the one-sided filtered Fourier sums. 

The same filter is also used to other functions such as u(x) = e co *( 1-5x ), u(x ) = 
etc., with similar results. We would also like to point out that even if we have restricted 
our discussion to the Galerkin method (1.2), the collocation case (1.3) can be analysed in a 
similar fashion with some assumptions on the location of the discontinuity between the grids 
(say, in the middle of two grids). We have performed the same numerical experiments with 
collocation method, and have obtained similar results. Details are omitted here. 

The main potential of one-sided filters is in solving hyperbolic partial differential equa- 
tions with discontinuous solutions. They can either be used during the final stage of post- 
processing to recover accurate point values of the numerical solution near the discontinuity, 
or used in each time step to obtain non-oscillatory numerical solutions. In order to per- 
form the latter one might obtain three Fourier sums, namely right-sided filtered u^} R (x), 
two-sided filtered u^} T (x) and left-sided filtered u^ L (x), at each point (this only involves two 
additional FFT’s if implemented in the phase space (1.5)), and use some local criterion to 
decide whether there is a discontinuity nearby, and if yes, to which side, then to decide which 
filtered solution to use. This is very similar to the ENO idea in finite difference [5], and is 
currently under investigation. 
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Figure 1: The sawtooth function (background solid line); the Fourier sum (1.1) with N = 16 
(solid line) and the two-sided filtered Fourier sum using (1.8) with p = 6 (dashed line). We 
can see that the two-sided filter removes oscillations away from the discontinuity but still 
leaves over- and under-shoots near the discontinuity. 



Figure 2: The right-sided filters : the real parts Re(a £ ) (solid lines) and the imaginary 
parts /m(of ) (dashed lines), for (a) N = 8; (b) N = 16 and (c) N = 32. 


Figure 2(a) 
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Figure 2(b) 



Figure 2(c) 
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Figure 3. The right-sided filters cr^\ the magnitudes (solid lines) and the arguments 
Ar 9(*k) (dashed lines), for (a) N = 8; (b) N = 16 and (c) N = 32. 

Figure 3(a) 



Figure 3(b) 
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Figure 5: The function (3.4) shifted by 7 r (solid line) and the one-sided filtered Fourier sum 
(the plus signs), with N = 16. 



Figure 6: The errors |ti(x) — u^(a:)|, of the one-sided filtered Fourier sum to (3.4), in logrithm 
scales, for N = 8, 16, 32. 
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